Semiclassical framework for the calculation of transport anisotropies 
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We present a procedure for finding the exact solution to the linear-response Boltzmann equa- 
tion for two-dimensional anisotropic systems and demonstrate it on examples of non-crystalline 
anisotropic magnetoresistance in a system with spin-orbit interaction. We show that two decou- 
pled integral equations must be solved in order to find the non-equilibrium distribution function 
up to linear order in the applied electric field. The examples are all based on the Rashba sys- 
tem with charged magnetic scatterers, a system where the non- equilibrium distribution function 
and anisotropic magnetoresistance can be evaluated analytically. Exact results are compared to 
earlier widely-used approximative approaches. We find circumstances under which approximative 
approaches may become unreliable even on a qualitative level. 

PACS numbers: 03.65.Sq, 85.75.-d, 75.30.Hx 



I. INTRODUCTION 

The change of electric resistance upon varying 
magnetization direction is an old and well-known 
phenomenor^ii^i^ with applications in spintronics.'*'^ Al- 
though the experimental observation of this anisotropic 
magnetoresistance (AMR) is rather direct — first ac- 
complished as early as 1857 — its theoretical under- 
standing is far from being complete. It has long been 
clear that the phenomenon arises from the combined ef- 
fects of magnetization and spin-orbit interaction. Dis- 
regarding the crystalline anisotropic background, the 
magnetization-broken symmetry between two chosen di- 
rections and unequal resistivities along these has been 
described within different models. In transition metal 
ferromagnets, the anisotropy was ascribed to asymmetric 
scattering due to different parts of the spin-orbit interac- 
tion Z-S* = ^{L^S+ + L^S-)+LzSz and the mechanisms 
considered were dubbed the L-S+ model, ^ LzSz model, ^ 
or a combination of both.'^'S Later, when computational 
power became sufficient for such task, ab initio calcula- 
tions were performed^ for permalloy and reached a good 
agreement with experiments. However, no direct link 
between the ab initio and the model calculations listed 
above has been established, probably due to rather com- 
plex band structures involved. On the other hand, such 
link between microscopic calculations and a simple phys- 
ical model was recently found in the diluted magnetic 
semiconductor-iSiii (Ga,Mn)As whose band structure is 
much simpler. 

Despite the long history of the AMR research, the 
question has not been answered to date of how a rigor- 
ous transport formalism for anisotropic systems should 
be formulated. Instead, the transport anisotropy has of- 
ten been discussed only in terms of the asymmetry in 
scattering amplitudes between two states on the Fermi 
surface. Current availability of materials with relatively 
simple band structure motivates the quest for more pre- 
cise theories of AMR. The present article discusses how 



the semiclassical Boltzmann equation should be solved 
in anisotropic systems, using an example of the model 
two-dimensional (2D) electron system. This allows us to 
put the previous approximations on rigorous grounds. 

The conductivity of a given solid in the regime of linear 
response to the electric field £ can be evaluated within 
the semiclassical picture once we have found the distri- 
bution function satisfying the Boltzmann equation. In 
the literature, this non-equilibrium distribution function 
is approximated in several ways. The relaxation time 
approximation (RTA) relies on calculating the transport 
relaxation time r from the scattering amplitudes w(fc, k') 
between two states on the Fermi surface using 
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w[k, k') [l — cosi^jT^:,] 



(1) 



For isotropic systems, where w depends only on the angle 
■d^j:, between k and k' , the integral ([J) does not depend 

on the direction of k and the RTA provides in fact the 
exact solution to the Boltzmann equationi^. The scatter- 
ing rate 1 /r depends only on energy and it is constant on 
the whole Fermi surface once the Fermi energy is fixed. 

For anisotropic systems, keeping Eq. ([T]) in use pro- 
duces 1/r that depends on the direction of k. The 
non-equilibrium distribution function constructed utilis- 
ing the RTA and Eq. ([T]) may capture some aspects of 
the transport anisotropies but it is certainly not precise. 
This approximative approach underlies for example our 
previous calculations^^ or those of McGuire and Potter— 
and we refer to it as to the "l/r approach". 

An improvement was proposed by Schliemann and 
Loss.^'"* In what we will call the "l/r" & l/r^ approach" , 
they use, according to their notation, Eq. ^ to calcu- 
late I/t'I (fc), and they provide an explicit formula for the 
non-equilibrium distribution in terms of this quantity and 
of 
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In our article, we argue that in a general case the non- 
equilibrium distribution function cannot be exactly cal- 
culated by just evaluating two integrals such as Eqs. (|l|2p 
for each fc-point of the Fermi surface separately. In- 
stead, an integral equation must be solved that deter- 
mines the whole non-equilibrium distribution at once. In 
Section [ni we describe this exact "integral equation ap- 
proach" to transport in anisotropic 2D systems and then, 
m Section [ml we use a simple model system to illustrate 
how the procedure works. For this purpose we intro- 
duce the Rashba Hamiltonian combined with a scatter- 
ing potential due to randomly distributed charged and 
ferromagnetically ordered impurities. In this model, the 
AMR results from the spin-orbit coupled band structure 
and the broken time-reversal symmetry of the scattering 
potential. ^"'^ We explicitly calculate exact non-equilibrium 
distribution functions for several specific realizations of 
this model, starting from the ones with simple solutions 
and then proceeding to the more complex case. Through- 
out Section IIIII we compare our distribution functions 
and AMRs to results of the other two approximative ap- 
proaches. Section IIVI concludes the main body of the 
article by discussing the relevance of our model calcula- 
tions for the AMR in magnetic semiconductors and by 
summarizing the key elements of the theoretical frame- 
work we have developed. The Appendices contain details 
of our calculations and also a more thorough description 
of the 1/t and 1/rll & l/r-"- approaches. 



II. THE FRAMEWORK 

Our central goal is to obtain the distribution function 
/ = /(fc, £) of a conductor displaced from equilibrium by 
a small bias represented by a weak homogeneous electric 
field £. We start with the familiar form of the Boltzmann 
equation in 2D 



e£-v{k)[-d.h) 



^nj{Kk')[f{k,£)^f{k',£)] 

(3) 

for a steady state of a spatially homogeneous system. 
Derivation of this equation is described for instance in 
Ref. Equation ^ is valid up to linear order in 

\£\ and it assumes small deviations of f{k,E) from 
the equilibrium distribution /o — fo{k). The velocity 
V = (l//i)Vfce^ is implied by the band dispersion e^, and 
e is the charge of carriers. The scattering rate w (per unit 
area of the reciprocal space) from the state k to k' needs 
to be specified according to the microscopic origin of the 
scattering; specific examples can be found in Appendix HSI 
or in Eq. (38) of Ref. [ij- For now we only assume that 
the scattering is elastic, that is w{k,k') oc S(^ej: — e^,). 
Focusing on the AMR, we do not include anomalous 
termsi^ii^ like the coordinate shift related to the side 
jump in the anomalous Hall effect into the right hand 
side of Eq. ([3]). Equation ([3]) can be applied to multi- 



band systems where k is replaced by a compound index 
containing the wavevector and band index and the inte- 
gral by integration over the wavevector and summation 
over bands. 

The solution to Eq. ([3]) is a function both of k and 
£. Focusing first on the latter, we can write it as Taylor 
series 



f{k,£) = h + £,deJ + £ydeJ- 



E 



£i.£jd£^d£.f 



(4) 

Being interested only in the linear order of the electric 
field components {£x, £y), we will truncate the series after 
the third term. For the simplicity of notation, we will 
now assume that the band structure is isotropic, ej: — Ck 
and —Sefo — S{e — e^). The anisotropy can still pervade 
into the transport via w, for instance due to anisotropic 
scatterers. The more general Boltzmann equation for 
anisotropic bands is treated in Appendix I Gl 

We define two angles 4>, 9 as £ — £{cos9,sm9) and 
k — A:(cos (/), sin 0) and Eq. ([4]) becomes 

/(0, 0)-fo= £{A{<I>) cos 9 + B{(j>) sin 9) (5) 

after the truncation, where A{(f)) = dg^f and B{(f)) — 
dg f. The non-equilibrium distribution is now expressed 
in terms of two functions of (f> which must, according to 
Eq. ^ with Eq. (O inserted, fulfil 



(6) 



w{(t>)a{(l))- j d4)' w{(j),<j>')a{(l)' 



sm ( 



We define here ^(0) = -ev{-d^fo) a{(l)), B{(j)) = 
—ev{—defo)b{(j)), and w{(j)) = J del)' w{<j),(j)'), where 
w((/), (/)') — (27r)^^ j k' dk'w{k, k') now includes the origi- 
nal transport scattering rate w{k, k') and also the density 
of states. 

The integral equation ^ with two variables (j), 9 can 
be decomposed into two independent integral equations 



coscf) = w{(j))a{(t)) - J d(f>'w{(j>,(t)')a{(j)') (7) 
sincj) = w{(j))b{(l))~ I d(f>'w{(j),(l)')b{(j)') (8) 



whose solutions a((/)), b{(f>) inserted into Eq. yield the 
exact solution of the Boltzmann equation ([3]) up to the 
linear order in £. 

The two decoupled inhomogeneous Fredholm 
equations^^ of the second kind (|7l8p can be most 
conveniently solved in terms of Fourier series. For 
special choices of wlcj), cj)') or band structure anisotropy 
(see Appendix l G[) . the series may contain only few terms 
and reduce to an ansatz for /(0, 9) such as Eq. (15) in 
Ref. 17. We explain the general procedure how to solve 
Eqs. (|7l8p on three examples below. 
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III. THREE EXAMPLES WITH RASHBA 
SYSTEM 

To illustrate how the above outlined procedure works, 
we choose the 2x2 Rashba Hamiltonian^ in the basis of 
plane waves 



(9) 



where ax.y are the Pauli matrices, and A is the Rashba 
parameter. In addition to Eq. we consider scattering 
on dilute charged magnetic impuritie o^'^'^^ described by 
the operator V, 



V/Vo = a + cTx = 



a 1 
1 a 



(10) 



that is impurities containing short range electric and fer- 
romagnetically ordered magnetic potentials. The quan- 
tity a is the (dimensionless) strength of the electric part, 
relative to the magnetic part, of the 'electro-magnetic 
scatterer' whose magnetic moment was chosen to be 
along the x direction. The magnitude Vq and other as- 
pects of this model are discussed in Section IIVI and Ap- 
pendices El and [B] 

We now calculate the non-equilibrium Boltzmann dis- 
tribution function f{k, £) = f{(j), 9) for this model in sev- 
eral special cases. To facilitate relevant comparison be- 
tween the 1/r and l/r" fel/r-*- approaches and the ex- 
act integral equation approach of Section [ill we calculate 
/((/), 6) and evaluate the AMR within the approximative 
approaches as well. 



A. Single band and magnetic scatterers 

The first special case of the model above concerns 
purely magnetic scatterers (a = 0) in the situation 
when the Fermi energy cuts the spectrum of the Rashba 
Hamiltonian ([9]) precisely at the A: = degeneracy point 
(ef = in Fig. [T]). We further disregard this single point 
of the Fermi surface and consider only the '+' band. This 
case offers the simplest way to explain the calculation of 
/ outlined in Section HI! 

The dimensionless scattering probability correspond- 
ing to V of Eq. (dni) is 



w{c^,cj)')/K=-[l 



5(0 + 0')] 



(11) 



This result, including the dimensionful prefactor K, is de- 
rived using the Fermi golden rule in Appendix 1"^ and 1^ 
Although w{4),(t>') does not explicitly depend on the 
Rashba parameter A, the presence of the spin-orbit cou- 
pling, combined with the symmetry breaking scattering 
potential, has the crucial implication that w depends on 
absolute values of angles and 0'. This leads to the non- 
zero anisotropy of the magnetotransport, in contrast to 
the isotropic case in which w{(j)^ <f>') depends only on the 




FIG. 1: Summary of the Rashba model defined by the Hamil- 
tonian ((9|. The three-dimensional plot shows energy disper- 
sions; its cross section along one (arbitrary) direction in the 
fc-space is also shown. The spin textures for the pair of eigen- 
states and |fc— ) for each k at the respective Fermi sur- 
faces with £f > are shown on the top. Unlike the diameters 
of the -|- and — Fermi surfaces, the eigenstates (spin textures) 
are independent of the Fermi energy ep, according to Eq. (|42|l . 



relative angle <j) — 4>' between the incoming and outgoing 
momenta. The total scattering probability, implied by 
Eq. nil), reads 



w(0) = Kit . 



(12) 



Note that despite the independence of w on in the 
special case considered in this subsection the resulting 
relaxation times and conductivity are indeed anisotropic. 

We will now look for the solution a(0) to Eq. ([7]) in 
the form of Fourier series 



a(0) 



Oci COS ( 

flsi sin (, 



' ac2 COS 2(j) - 
sin 2(j> - 



(13) 



Owing to the trivial form of w{4i, cf)') (of its Fourier spec- 
trum) the integral in Eq. ([7]) can be readily calculated 
and Eq. ^ assumes the following form: 



1 



■ cos ( 



— flci COS (j) + ac2 COS 2(f) + cos 30 + . . . + 



-asi sm( 



sin 20 + 0^3 sin 30 - 



The only non-zero coefficients in the Fourier series (|T3)) 
are therefore oq and Uci- The solution of Eq. ([7]) then 
reads 



a(0) = ao ■ 



2/3 



(14) 



Conservation of the number of particles requires oq to be 
zero. 

A completely analogous procedure applied to Eq. ([8]) 
yields a system of equations for coefficients bg, bgi, bci 
which give 



6(0) 



ttK 



(15) 
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The complete solution up to linear order in £ to the 
Boltzmann equation written using Eq. ([5]) is there- 
fore 



f{<j),e)^fo~ev£{-d,fo) 



ttK 



1 



Let us now compare this result to the approximate 
approaches outlined in the Introduction. The non- 
equilibrium distribution in the l/r approach is (see Ap- 
pendix [C| 



/(0,0)-/o = -ev£{~dJo)— X 

ttK 



(17) 



cos ( 



sm ( 



cost 



3-2 sin^ 



smf 



3-2 sin^ 



while in the I/t" & approach (see Appendix rPjl. we 
obtain 



/(^,0)-/o = -ev£{-dJo)^ X 

ttK 



(18) 



COS0 



3 cost 



2sin^ 



COS( 



-|- sin t 



9 -h 4 sin'' 
,3sin(^(l 



S sm 



cos< 



2sin'^ 



4 sin" 



Distribution functions in Eqs. (|16I17I18|) are significantly 
different. To quantify the differences, we use these three 
distribution functions to calculate the AMR, defined as 



AMR = 



'yy 



(19) 



'yy 



and having the meaning of the (relative) difference in re- 
sistivity for current flowing parallel and perpendicular to 
the direction of the scatterer's magnetic moment, respec- 
tively. The conductivities are calculated from the current 
implied by the non-equilibrium distribution f{k,£) 



ev{k)f{k,£) 



(20) 

7r/2)/£ where 



i.e., as Uxx = j{d = 0)/£ and cjyy = j{ 
j{e)=Jdq>f{cl),e)\v\ cos(0-0). 

The AMR value of 1/2, obtained from the exact 
non-equilibrium distribution function in Eq. (|16|) . is 
markedly different from the results of the approxima- 
tive approaches. The l/r approach underestimates the 
AMR by almost a factor of two (AMR « 0.27), and 
the l/rll&l/r^ approach predicts even a wrong sign 
(AMR« -0.11). 

Before we proceed to comparing the three approaches 
on other realizations of our model disordered 2D system, 
let us make a remark about the distribution functions 
above. The non-equilibrium part of the distribution func- 
tion in Eq. ([TT]) was obtained as —ev-£{~defo)T{<j)) with 
v-£ = v£ cos{e - 0) and t((/.) = {2/TrK){3 - 2 sin^ (jj)-^ 
as derived in Appendix [C\ Analogous factorization of 



the bracket in Eq. ((T8|) or Eq. (|T6|) is not possible, re- 
flecting the fact that no scalar relaxation time can be 
attributed to a given fc-state in these approaches. How- 
ever, the I/t'I &1/t^ approach still unambiguously as- 
signs relaxation-rate-like quantities, a pair of (not neces- 
sarily positive) values 1/Tll(fc), 1/T^(fc), to each A;-state, 
independent of the electric field direction (determined by 
9; see Appendix l^j) . It remains an open question whether 
also the exact solution of the Boltzmann equation, such 
as Eq. pop , can be meaningfully interpreted in terms of 
6'-independent quantities related to scattering. 



B. Single band and electro-magnetic scatterers 

We now extend results of the previous section by relax- 
ing the condition a = 0, that is we consider the complete 
scatterer with electric and magnetic parts of its potential 
added up coherently, as defined by Eq. pU)) . The exten- 
sion is straightforward although the algebra involved is 
richer than for the previous model. The dimensionless 
scattering probability w((/), (t>')/K and w(0) are 

w{(j),(j)')/K = i[l-cos(<?!) + (/)')+"^(l + cos(?!)-0'))] 
+a(sin (j) + sin (f)') 

w{(j)) = TTK{l + a'^ + 2asm(t)), 

(21) 

as shown in Appendices [X] and [B] Note that w{(j)) is 
no longer constant, which is here the direct reason of 
the more complex algebra needed. We again look for 
the solution of Eq. ([7]) in the form of Fourier series 
and find that the higher order coefficients ac2, Ocs, . . . are 
now no longer zero. Instead of Eq. (|14p . we get a system 
of an infinite number of linear equations which is not 
surprising, given that Eq. ([7]) is an integral equation in 
its general form. 

This system of equations can be solved using a parti- 
tioning method, described in Appendix [Ej Herein, we 
segregate the variables into three groups: {oq, Cd, a^i}, 
{ac2, ac3, ■ • ■}, and {0^2, 0^3, . . .}. The first group must 
obey 



l/{TrK) = (1 + a'^)aci + 0:0,2 + ^(1 - a^)aci 
= {1 + a^)asi - aac2 - hdsii^ + a^) ■ 



(22) 



and flo = as in the previous Subsection. Equations (|22p 
originate from comparing the coefficients in front of the 
cos (j) and sin terms of Eq. ([7]) with Eq. inserted. 
Separate treatment of the other two infinite systems of 
equations, described in Appendix 1^ yields 



ac2 



as2 



—ttcia for \a\ < 1 
— Cci/a for |q;| > 1 



(23) 



Together, Eqs. (|22l23p form a closed system for Od, and 
flsi which thus read 



1 r4/(3-a2) for |a| < 1 
" 2^if \ 4/(1 + a2) for|a|>l 



Osl = . 
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The solution to Eq. ([7]) for |a| < 1 is then 
1 4 



2t:K 3-a' 



■ cos (f) + ac2 cos 20 + as2 sin 20 - 



(24) 

For the evaluation of current and AMR using Eq. (fT9|) 
there is no need to know the higher order terms of 9) 

(by virtue of J^^ cos cj) cos 2</) = etc.). However, keeping 
all higher order terms of Eq. (I13p in the derivation was 
necessary for obtaining the xscorrect form of Eq. (|^^ and 
also correct expressions for constants ad and a^i at the 
end. 

We again repeat the same procedure for Eq. ([8]) , obtain 
&((/)), and finally we complete the calculation by writing 
down the non-equilibrium distribution function: 



/(</>, ^)-/o 



-evSi-dJo)— X 
ttK 



cos 9- 



(25) 



sm ( 



3 - a2 ■ ' 1 

for \a\ < 1, while for |q!| > 1 the bracket is replaced by 



COS( 



sm ( 



cosf 



1 



1 



(26) 



The dots symbolize cos 20, sin 20 and higher order terms 
which as emphasised above do not contribute to the 
AMR. The divergence of this expression for jaj ^ 1 will 
be discussed in Section HVl 

Evaluating the AMR using the distribution func- 
tion (|25l26p and Eqs. (|20ll9p amounts to comparing the 
coefficients in front of the cos 6 cos and sin 9 sin sum- 
mands. We get 



AMR== 



1 



lal < 1 



AMR = 4t 
a'' 



lal > 1. 



(27) 

We conclude the study of the single-band model by 
comparing this AMR to the results of the approximate 
l/r and l/rll&l/r^ approaches shown in Fig. [UJa). 
While the l/r approach can be regarded as only quan- 
titatively inaccurate, as already suggested by the results 
of the previous Subsection, the apparently more sophisti- 
cated l/r" & I/t-^ approach yields remarkably large de- 
viations from the exact AMR. 



C. Two bands and electro-magnetic scatterers 

We now consider the case when the Fermi energy is 
above the fc = degeneracy point of the Rashba bands. 
Let us first explicitly write down the scheme of Sec.HIlfor 
a two-band system. Considering distribution functions of 
the and '— ' bands, denoted by /+ and /_, Eq. ^ is 
replaced by two coupled equations 



e£ ■ v+{k) — 



de 



(28) 
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FIG. 2: AMR as a function of the relative strength, a, of the 
electric and magnetic components of the impurity potential 
(see Eq. (|1Q|) for the definition of a). Dashed and dotted lines 
denote results of the approximative l/r and l/r" & l/r"'" ap- 
proaches, solid line shows the exact result of the integral equa- 
tion approach, (a) Single-band model, (b) two-band model 
with A ^ 0. 



[{w++{k,k')+w+^{k,k'))f+{k)- 
-w++{k, k')f+{k') - w+^ik, fc')/-(fc')] 



and 



-e£ ■ V- (fc) 



<Pk' 
(2^ 



[(«;_^(fc,fc')+^-+(fc,fc'))/-(fc)- 
-w..(k,k')f.{k')-w.+ {k,k')f+{k')\ 



The scattering rate w ji now also bears the indices of the 
initial (z) and final (/) band. We will abbreviate the 
equilibrium distributions /o(ej_|_) by /o±. 

Assuming isotropic bands 6^+ and efe_, we seek a so- 
lution of Eqs. ((28)) in the form of 



/+( 



) - /o = -e£v+ {-dJo+){a+ (0) cos 9 + b+ (0) sin 9) 
) - /o = -e£v- i-djo- ) (a- (0) cos 9 + b- (0) sin 9) 



and the four functions a±{(j)), h±{(j}) must fulfil 



(29) 



u;+(0) a+(0) 

d0' [u;++(0,0')a+(0') ('/'')] 
iD_(0)a_(0) 

-p' [w — (0, 0')a_(0') w_+(0,0')a+(0')] 



6 



sin0 = w^{(f))b+{(j)) 

- J dc^' K+(0,0')&+(f) + 'i'+-('^,0')&-(0')] 

siii0 — w-{if)) b^{4>) 



(31) 



b' <j>')b^ (0') + w^+{cf>, <j>')b+ (0')] 



where ix;/i(0, (/)') = (27r) J k' dk'w fi{k, k') a.nd ■Wi{(p) — 
J d(|3'[w+^{(t),(t)') + w-i{(t), (!>')]. Note that Eqs. §0^ are 
decoupled from Eqs. ([51]) . 



J. Evaluation of f+ and f- 

The dimensionless scattering probabihties for the 
complete electro-magnetic scattering operator given by 
Eq. Uni) are 



W{++/—)i(f>',(f>)/K = 
= i [l - cos{(f) 



COS 



'1 



(+/-)Q(sin0 - 



= i [1 + cos(?!) + + - cos(0 - (/)'))] + 

-|-Q;(sin (j) — sin 0') 
(32) 



and w±{(f>) = 2ttK{1 + q ± 2a sin ( 



For simplic- 



ity, we assume that the constant K (and the density 
of states, as explained in Appendix fG]) is the same for 
both bands. This occurs in the Rashba model ^ in 
the limit oi Xkp ^ cf and we will call this the A — > 
limit. Details of the derivation of Eq. ([32| are given in 
Appendix [A] and [B) 

In a close analogy to the single-band case, equa- 
tions ([SO]) produce two coupled infinite sets of linear equa- 
tions for coefficients of 



0'±i4') — aQ± + aci±cos(t) + ac2±cos24> - 
+ asi± sin (j) + as2± sin 2(p + 



(33) 



These may again be reduced to two coupled 3x3 sys- 
tems for variables ao±, Ocii, Osii using the partitioning 
method. Mathematically, their solution 



^0+ ^ flo- 

0-C1+ = flcl- 

0-31+ = flsl- 



= 

_ 1 

a 

= 



2ttK 



a 

1/a 



for |a| < 1 
for \a\ > 1 



(34) 



leaves ao+ -I- ao- undetermined and, physically, particle 
number conservation again dictates that this constant is 
zero. Terms in Eq. ((33| containing higher multiples of 
(j) are again not contributing to the current and to the 
AMR but their coefficients can be evaluated within the 
partitioning procedure. 

Applying the same procedure to Eqs. (f3T|) leads to 

^±(0) = 6o± + 6cl± COS (/) + 6c2± COS 20 + . . . -I- 
-|- bsi± sin (j) + bs2± sin 2(1) + . . . 



with 










—Za 




-bo- 


1 + a2 


bci+ 


bci- 


= 








bsi+ = 


bsi- 


= bs = — 



- 1 2ttK 



(-1) for \a\ < 1 
for \a\ > 1 



The non-zero value of bo+ — bo- means that the scattering 
redistributes particles between the two bands. In another 
system, where the two bands would have different net 
spin polarization, such redistribution would correspond 
to the polarization of the particles by impurities. The 
overall particle number conservation nevertheless again 
requires 6o+ + ^o- = 0. 

The two non-equilibrium distribution functions are 
now 



/±( 



/o - e£v{-dJo±) 



1 



2ttK 



(35) 



1 



cos y COS ( 



1 



■ sm f sm ( 



± 



■ sm( 



for \a\ < 1. The distribution function for \a\ > 1 is 
given by Eq. (j35p with the term in the square brackets 
multiplied by l/a^. 



2. Comparison to the approximate approaches 

The explicit calculation outlined in Appendix FfI shows 
that again all coefficients appearing in front of the co- 
sine terms of Eq. p3p are non-zero. The infinite series, 
however, can be summed up and the complete exact non- 
equilibrium distributions fulfilling Eqs. (|28p read 



f±{(b,e) = h-e£v{-d,h±) 



1 



2ttK 



(36) 



,1 



±a + sin ( 



1 + ± 2a sin c 



1 — a^ 1 + a^ ± 2a sin ( 
. „ ±a ' 



smt 



1 



a^ 



for |a| < 1. The 1/r approach (see Appendix fC)) leads 
to a similar but not identical approximative result 



/±(0,^) = k-e£v{-dJo±) 



1 



2t:K 



(37) 



COS( 



sm ( 



COS! 



1 + a^ ± 2a sin ( 



1 + a^ ± 2a sin c 



and the l/rH &1/t-'- approach gives precisely the same 
result as Eq. p7|) because vanishes in the two-band 
case (see Eq. (|50l) in Appendix [D| . 

Remarkably, the AMR calculated from the distribution 
functions of Eqs. or and of Eq. i.e. the 
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exact and the two approximative results, comes out to 
be the same 

AMR = , |a| < 1 AMR = 4r , |a| > 1 . (38) 

Plot of this function is shown in Fig.[2l|b). 

For the two-band Rashba model with small A, we thus 
conclude that the discrepancy between the exact and ap- 
proximative approaches remains only on the level of the 
complete non-equilibrium distributions (in the higher or- 
der terms that do not contribute to the current). We 
speculate that the equal results for AMR were not ob- 
tained by coincidence but because the A — + system has 
a higher symmetry than the single-band model which is 
chiral. These symmetries are briefly commented in Ap- 
pendix [d] 

IV. DISCUSSION AND CONCLUSION 

Let us start this discussion section with a remark on 
results shown in Fig.[5{a). The AMR takes on a singular 
value of 1 at a = 1 in all three approaches. This reflects 
the (1 — a^)~^ divergences of all non-equilibrium distri- 
bution functions (see Eq. (^5)1 for example). The origin 
of this divergence is as follows: the scattering operator V 
in Eq. pop with a — \ can annihilate one particular state 
|fc-|-) on the Fermi surface, as seen from Eq. (|^^ and the 
spin textures in Fig. [TJ The state has its spin aligned par- 
allel to the moment of the magnetic impurities, i.e., along 
the i-axis. For the Rashba model this implies that the 
/c-vector of this state is parallel to the y-axis, more pre- 
cisely 4> = — 7r/2. Within the l/r approach, the fact that 
(a -|- cra;)|A:-|-) = then implies that this state has an infi- 
nite (transport) relaxation time as dictated by Eqs. ((40)) 
and ([31]) of the Appendix. Consequent calculation focus- 
ing also on other states |fc) contributing to the current 
shows that this singularity is strong enough to produce 
(jyy oo when \a\ 1. 

The current (and AMR) calculated for |a| close to 1 
are clearly inconsistent with the linear-response basis of 
our theory approach (incorporated in Eq. ([3])) and are 
therefore not physically relevant. On the other hand, the 
impurity operator pU]) is idealised compared to realistic 
systems where the electric and magnetic part of V will 
depend at least slightly differently on fc. This modifica- 
tion suffices to remove the singularity in conductivity. 

Pointing our attention more towards experiments, let 
us now discuss the relevance of the Rashba model with 
dilute charged magnetic scatterers. Our original moti- 
vation comes from the study of the diluted magnetic 
semiconductor^^ (Ga,Mn)As. Mn atoms, when substi- 
tuting for Ga, introduce both the magnetic moments 
and holes to the material. The former via its d-electrons 
and the latter because their valence number is one less 
than that of Ga. An 'electro-magnetic scatterer' model 
as defined by Eq. ([TU)) is therefore relevant to describe 



the Mn atoms which constitute by far the most frequent 
source of scattering in (Ga,Mn)As. Indeed, it is possible 
to qualitatively explain trends for AMR in (Ga,Mn)As 
based on this model of scattering and by neglecting the 
exchange splitting of the (Ga,Mn)As spin-orbit coupled 
valence bandj ^°i^^ 

The Rashba model employed in this article provides ar- 
guably the simplest unpolarized spin-orbit coupled band 
structure in which the anisotropic scatterer mechanism 
fully determines the AMR. The other two mechanisms, 
which are the anisotropy of the group velocity and the 
anisotropy of wavefunctions of the spin-split spin-orbit- 
coupled valence band and which only quantitatively mod- 
ify the calculated AMR in (Ga,Mn)As, are completely 
absent in this model. The simplicity of the present 
model relies mostly in that it is 2D and it considers 
two rotationally-symmetric rather than six warped bands 
of (Ga,Mn)As. The integral equation approach can be 
straightforwardly extended to (Ga,Mn) As or other three- 
dimensional systems with more (n > 2) bands. However, 
the calculational complexity will be considerably higher; 
the two functions a±((?!)) of one variable will be replaced 
by n functions of two variables (two angles parametrizing 
the Fermi surface in three dimensions). 

Turning attention towards possible experiments, the 
calculations presented in this article are most relevant 
to asymmetric n-type heterostructures doped with mag- 
netic donorSfi^i By changing the Fermi level via doping, 
the effective strength a of the electric part of the scatterer 
should change because the scattering amplitudes depend 
on the Fermi wavevector which is a typical measure for 
involved momentum transfersi^ Consequently, by polar- 
izing the magnetic moments in-plane, the AMR defined 
in Eq. p9)) should be measurable and follow predictions 
shown in Fig. (Ha). 

An alternative to doping by magnetic donors is to use 
an n-type heterostructure co-doped with magnetic im- 
purities. Experimental study of a IILV or II- VI het- 
erostructure with dilute Mn doping and heavy remote 
n-doping could be revealing. Depending on the magnetic 
impurity character (either acceptor or neutral), by vary- 
ing the Fermi level, we could again effectively change a 
and/or interpolate between the single-band case (tp = 0) 
and the two-band case [tp ^ \kp). The challenge in this 
experiment would be to keep the scattering on Mn the 
dominant (or at least strong) mechanism of relaxation. 

Rather than these experimental suggestions, however, 
the main message of this paper should be of theoretical 
character. We have presented a framework to calculate 
exactly the conductivity in anisotropic systems within 
the semiclassical linear-response theory. This procedure 
was demonstrated on three simple and analytically solv- 
able models. We found that in some special cases of 
high symmetry the previously employed approximate ap- 
proaches may yield the same AMR as our exact theory. 
In general, however, only the exact non-equilibrium so- 
lution to Boltzmann equation of the form of an integral 
equation over the whole Fermi surface, rather than of 
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effective scattering rates at each individual fc-point in- 
dividually, provides a reliable account of the anisotropic 
transport. 
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APPENDIX 

A. Scattering rates 

We evaluate the scattering rates using the Fermi golden 
rule. Probability Wfi of transition between states \i) 
and I/), induced by a perturbation described by time- 
independent operator V, equals 

Wf.^^\{f\n)\'S{ef-e,), (39) 

where ej/i is the energy of the final/initial state. 

Considering many scatterers described by the operator 
V distributed randomly with areal density n.^, the scat- 
tering rate per unit reciprocal space between the k- and 
fc'-state equals 

wik,k') = ^n,V,'\(k'\V/Vo\k)\'S{ej:~ej:,) (40) 

within the lowest order of the Born approximation; con- 
trary to the case of the anomalous Hall effec t "'^^1^° , this 
order of the Born approximation is sufficient for the cal- 
culation of the AMR. Note that the dimension of the 
scatterer strength Vq is Jm^, making V /Vo dimension- 
less in the Fourier space. 

Finally, assuming isotropic parabolic dispersion — 
h^k'^/2m,, the density of states equals ■m/{hh) per spin, 
so that 



w{(l),(l)') = 



1 



k' dk'w{k,k') = (41) 



:\{k'\V/Vo\k)\'^K\{k'\V/Vom' 



This is the definition of the dimensionful constant K used 
in Eq. (fTTj) and later on. Its value determines the absolute 
value of conductivity but it cancels out in the definition 
of the AMR, see Eq. ^ 



B. Scattering matrix elements 

We calculate the matrix elements of the scattering op- 
erator V in Eq. (fTU]) with respect to the basis 



V2 ( T«e** ) Vl'' 



k-r 



(42) 



where k — fc(cos 0, sin (p) and A is the system area. Vec- 
tors and |fc— ) are the eigenstates of Hamiltonian ^ 
with eigenvalues ti^ k"^ / 2m— X\k\ and fi^k"^ /2m+X\k\; their 
(expectation value of) spin a = (axjCTy) is illustrated in 
Fig. [TJ The scattering operator V in Eq. ^TU\i is expressed 
in the basis of plane waves times spin up and spin down 
states. It does not depend on fc, k' so that it corresponds 
to short-range impurities ((5-scatterers). For a — > cx), this 
would be a non-magnetic charged impurity of strength 
aVb, and for a = it is a purely magnetic impurity of 
strength Vq. 

Owing to the (5-scatterer character of V, the matrix 
elements of V/Vq in the basis ((42|) depend on k only 
through (f) and not through k = \k\. We take k = k' and 
get 



1 



(fc'+|a + cr^|fc+) 

(^'-|a + (7^|^-) = i[«e"^-2e-"^' +a(l-he'('^-*'))] 

(A?'+|a + (7^|fc-) = i[ie*'^ + ie-"^' +a(l -e'^-^^"*'))] 

= (fc-|a + CT^|A?'-h). 

Taking the absolute values squared leads using Eq. (|4T|) 
to Eq. ([311), to Eq. dH]) (the '++' element), and to 
Eq. ini) ('++' element with a = 0). 



C. The 1/r approach 

Non-equilibrium distribution function in an isotropic 
(eg = Cfc, and isotropic scatterer) two-band system can 
be shown to be 



/±(0,e) = /±( 



= fo-ev±-£ - 



a/, 



de 



T±, (43) 



where the relaxation times for -I- and — bands may de- 
pend on k only through energy e^. This fact, that for 
fixed energy the relaxation time as defined in Eq. ([T]) is 
constant, is a direct consequence of the scatterer isotropy 
w{k,k') = w{'dj^i^,). For clarity, we stress that in an ri- 
band system there are the total of 71^ scattering rates 
between pairs of bands, 



1 

Ua 



(Pk' 



;UH,a{k, k') 



\vb{k')\ 

l«a(fc)| 



COS "iJfTfT/ 



, (44) 



9 



that combine into n scattering times Ta, one for each 
band, according to the Matthiessen's rulsi^ 



Ta ^ Tba 



(45) 



We note that d^ij' measures the angle between Vb{k') and 
Vaik) but given the isotropy of the band structure, v{k) 
and k are parallel so that d^jij' = ^jif:,- Equation ([1]) is 
a single-band variant of Eq. (|44p for isotropic systems 
where v drops out. 

In the 1/t approach, we simply evaluate Eq. (|44)l for 
'w{k,k') ^ ^(''^fcfc/) and obtain A;-dependent 1/t. This is 
then inserted into the distribution function (|43p . losing 
thereby its property /±((/), 9) ^ f±{(j)-9). 

For the Rashba model with A ^ we get using Eq. (PT|) 
and Appendix FB] the following 



'(++/- 
''(+-/-+) 



^ ii'7r((+/-)2asin0 + 3- 2sin20 + a2) 
(0) =ii:7r((-/+)6asin(/>+l + 2sin2(/) + 3a2) . 

(46) 

The scattering rates in the two-band model are thus as 
simple as 



1 



1 



1 



l/(47ri^) 



'''±{<t') T+±{<t') T^±i4') l±2Q;sin( 



(47) 



This result, plugged into Eq. (^5]) . produces the non- 
equilibrium distribution function ([37]) within the 1 /t ap- 
proach. 

The relaxation time for the single band model is simply 
1/t(0) = 1/t++(^). Setting here a = leads via Eq. 
to Eq. p7|). 



better chirality), i.e. scatter more clockwise than coun- 
terclockwise — such as it is the case with skew scattering 
in the anomalous Hall effect, (iii) Contrary to the state- 
ment of Ref. the distribution function is not the 
exact solution to Eq. ((28)) for a general anisotropic sys- 
tem. The derivation of Eq. (|48|) presented in Ref. 13 is 
only valid if the expressions 



given by Eq. (27,28) of that reference are constant for 
each band (i.e. (/)-indcpcndent in our case). The most 
general distribution function this approach can therefore 
correctly capture must have the form 



/(0,0)-/o = Ci cos( 



C2 sin(0 - 61) 



while as the examples in Section |TT] show, the non- 
equilibrium distribution can have finer details than those 
of period 2tt in the angular variable (f> (and these details, 
when completely neglected, may even lead to wrong val- 
ues of the constants Ci , C2 above) . This original neglect 
of Ref. [H was later corrected by one of its authorai^ 
in the context of the specific Hamiltonian considered.— 
However, a general procedure for exact solution of the 
Boltzmann equation was not given. 

In our specific model, as described by the scattering 
matrix elements of Appendix [B] we get 



'(++/- 



L'(+-/-+)(0)] 

For the two-band model. 



Kn cos 0((-f /- 
Kn cos 0((-f /- 



)a - 



sin( 
sin ( 



(49) 



D. 



The l/r"&l/r^ approach 



1 



1 



1 



0, 



rf(0) 



0, (50) 



The prescription for the non-equilibrium distribution 
function suggested by Schliemann and Loss^';? can be sum- 
marized as follows: (a) evaluate the 'standard' formu- 
lae (|45l44p and denote the result as I/tJ ((/>); (b) cal- 
culate 1/t^(0) using formulae identical to Eqs. (|44l45p 
save the replacement of the bracket in Eq. (|44p by 
|t'b(fc')|/|i'a(fc)| sini9^^/); (c) write down the distribution 
function as 



f±i 



cos 



/o - e\v±\ \£\{-dJo±) X 



(48) 



-sin((/)-6l) 



T-L_ll 



Several remarks are in order, (i) Whenever van- 
ishes, Eq. (|48| simplifies to Eq. ([43]) of the 1/t approach, 
(ii) This 1/rll & 1/t^ approach is suitable for the descrip- 
tion of isotropic scatterers (the amplitude depends only 
on the angle between k and k' , the incoming and outgo- 
ing wave) which however may exhibit an asymmetry (or 



so that the l/r" & l/r^ approach reduces to the 1/r ap- 
proach in line with the comment after Eq. (pS)) . The 
single-band case, however, has a finite = T^_^_ so that 
the two approaches give different results. This is not sur- 
prising, since each Rashba band has a chiral spin texture 
but both of them together form a non-chiral pair, pro- 
vided they have both the same Fermi k (as it happens 
for A ^ 0), see Fig. [TJ The asymmetry of scattering 
expressed by 1 /r^ thus vanishes in our two-band model. 
To obtain the non-equilibrium distribution (jl8p within 



the l/r"&l/r^ 
Eq. (glD, 1/t^+ 
and expand cos( 
sinf?. 



.11 



of 



approach, we have to take l/r 
of Eq. ([iS]) . insert them into Eq, 
i — 9), sm(6 — 9) in terms of cos 9 and 



E. Partitioning method 

The actual infinite system of linear equations for 
variables oq, Od, flgi, . . . appropriate for the single-band 
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model assumes a structure suitable for partitioning if we 
perform the coordinate transformation ip — tt/2 — (j). We 
will now solve the integral equation Q using this coor- 
dinate (and use ip throughout Appendices 1^ and FF]) and 
transform the result back before we use it in Eq. (|22|) . 
The equation to be solved is now 



w{ip)d{ip) — / d(f'w{(f , (f')a{(f') = sin (f (51) 



with 



-|-a(cosi^ + cos 1^') 
iv{'p) = 7ri^(l + a'^ + 2q; cos 1^) . 

Inserting 

a(ip) =00+ del cos (p + . . . + dsi sin (p + . . . 

into Eq. ([5T|) . and comparing the coefficients at the con- 
stant, cosip, sini^, cos2ip, cos3ip, . . ., sm2ip, ... terms, 
we obtain the following infinite system of linear equations 
for do , Oci , Osi , , ac3 , ■ • • , as2 , as3 , • ■ ■ (in this order) : 

/ 















... 








... 










a 





... 








... 




* 


* 


* 








... 


a 





... 







a 





l + a^ 


a 


... 








... 














a 1 


+ 


a 








... 

















a 1 










... 











a 








... 


l + a2 


OL 


... 




















... 


a 



























: / 



V 

(52) 

The double line separates the left and right-hand side of 
the equations. The twelve asterisks in the first three lines 
of the system (|52|) correspond to the 3x3 system ([22l) , 
and the value of these coefficients will be unimportant 
within this Appendix. 

It is apparent that the system (|52p is almost block- 
diagonal. The partitioning method takes advantage 
of this structure and aims at solving three indepen- 
dent systems corresponding to groups {ao, ftci, Osi}, 
\dc2, dci, . . .}, and {as2, as3, . . .} of the original variables. 
The basic idea is to treat the only non-zero element of 
the off-diagonal block as a right-hand-side term. In ex- 
plicite terms, we rewrite for example the fourth and fifth 
equations of the system (|52|) 



\adc\ + (1 + a^)dc2 + aacs] cos2i^ 
\adc2 + (1 -f c?^dcz + adc/^ cos iCp 



as 



(1 -I- o?)dc-2 + adcz 

adci + (1 + o?)dc-i + aac4 



A 



The system of all 'cosine-term' equations of the sys- 
tem ([5^ (starting with cos2(^) can now be solved as a 
function of A. In other words, we are treating the central 
block of the matrix ([5^ . The still-infinite system to be 
solved is 



/l 



V 



a a 
a 





a 






a 



a\ 






(53) 



For the purposes of solving later the 3x3 system (f22|) . we 
in fact need to know only a part of the solution, namely 
ac2- To this end, linear algebra gives us a very quick 
answer. If we denote by D the determinant of the infinite 
matrix left from the double line in and by £)„ the 
determinant of the analogous n x n matrix, then if the 
system ((53)) were finite. 



ac2 = 



A_D^ 

Dn 



where the numerator equals the determinant of the nx n 
matrix left from the double line of ()53p with first column 
replaced by the column right from the double line. Con- 
sidering n — > oo, we immediatelly (after transformation 
ac2 — —dc2) get dc2 — —A as given in the first line of 
Eq. (|23p . This answer is, however, not completely cor- 
rect. 

The caveat of this procedure is that we should have 
been careful about taking the limit D = lim„_,cc-Dn- 
It turns out that the limit is finite only for \a\ < 1 
and then D = 1/(1 — a^) so that only in this case 
limDn-i/ Dn = (limZ3„_i)/(lim£'ji) which is obviously 
equal to one. The determinant D is infinite for \a\ > 1 
and only lim„^oo Dn-i/ Dn remains finite, namely equal 
to 1/a^ as one can readily see from the explicit formula 



Dn = l + 

In conclusion, we find 
ac2 



A for |a| < 1 
A/a^ for |a| > 1 



(54) 



and the transformation back from ip to cf) implies ac2 = 
— ac2 and dd = Gsi. 

Literally the same procedure works for the 'sine-term' 
equations of the system (|5^ . i.e. the lower-right block. 
The only difference is now that A = —ad^i and we use 
as2 = as2, dsi = ad. These two results, ac2,as2 with 
the corresponding definitions of A are summarized as 
Eq. (BSD. 

The key feature needed for this partitioning method is 
that A is a function only of dd (dd) and not of higher- 
order coefficients like 0^3. In this way, the 3x3 system 
of equations (|22p becomes closed after ac2 and have 
been inserted. 
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Finally, we stress, that if the coupling between the 
three subsystems had been neglected from the very be- 
ginning — this amounts to setting to zero the four el- 
ements in the off-diagonal blocks in the system ([5^ — 
the solution of the 3x3 subsystem represented by the 
asterisks would have been different. In this way, even 
though cos 2(j) and other higher terms do not contribute 
to the current calculated from the non-equilibrium dis- 
tribution their complete neglect from the beginning 
may produce wrong coefficients in the cos (f> and sin 
terms. 



F. Partitioning method — two bands 

In the case of two bands, we obtain two infinite systems 
of linear equations identical to the system ([5^ . one for 
variables with index, another for those with ' index, 
see Eq. ([55]) . Although the two systems are now coupled, 
the direct coupling exists only via variables ao±, Odi, 
dsi± corresponding to the upper left block. The parti- 
tioning method can therefore be independently carried 
out in the and '— ' sector. 

For all four infinite subsystems, of which the sys- 
tem (1551) is one, we obtain the almost the same result 



A(-a)"-2 for |a| < 1 
A/(-a)" for |a| > 1 , 



(55) 



for n > 2 and with appropriate definition of A for each 
subsystem, while acn-^ o,sn- obey Eq. (j55p with —a re- 
placed by a. All coefficients in the series ([55)) are thus 
non-zero. Nevertheless, Eq. ([33]) can still be summed up 
using 



E, . 1 + acos^ 
(—a) cos nip = 

n=0 



1 + a'^ + 2a cos ip 



and a similar formula for sines. We now transform back 
from if to 0, use A = =paaci± for the cosine ± parts of 
Eq. (I33p and A — ^adsi± for its sine parts, transform 
back dci± — flsii, asiifc = Ocii and finally get 



a+(0) = ao+ + 



C0S( 



J-C1+- 



+ asi+- 



a + sm ( 



I + a'^ + 2a sin (f> 1 + a'^ + 2a sin c 

cos 6 —a + sin 6 



Ocl- 



1 -t- — 2a sin ( 



1 -f — 2a sin < 



(56) 

Plugging the values of ao+, dci+i olsI+ from Eq. ()34p into 
Eq. (|56p . repeating an analogous procedure for the 6's in 
Eq. ([5T|) and inserting the results into Eq. ([^5)) , we arrive 
at Eq. ([55]) . 



G. Boltzmann equation in general 2D anisotropic 

systems 

Results of Section [nj were derived for a special class of 
2D systems where the band structure remains isotropic 
and the anisotropy is only introduced through the scat- 
terer and the scattering rate ?x;(fc, k'). 

The results of Eqs. (|7l8p for single-band or of 
Eqs. (|30l3ip for two-band system have to be slightly mod- 
ified for anisotropic 2D band structure. The wavevectors 
fc, k' of Eq. ^ or Eqs. ^ wiU stiU be bound to the 
Fermi level ep but their magnitude now depends on 0. 
That is, we have k — k((j)) and k = kfi = fc(cos0, sin0). 
We also tacitly assume that in each band and for each n 
there is only one solution k to ekn = cf- The calculation 
of w((f>,(l>'), compared to what is done in Appendix \~K[ 
becomes 



w 



(0,0') 



(27r)2 



k' dk'w{k,k') 



(57) 



-1 ^ 



The last expression should be understood as a function 
of (j>, (p' only; the derivative and k, k' are to be taken at 
the Fermi level, so that e.g. k' = fc'(cos 0', sint/)') and 

Further, the expression £ ■ v{k) in Eq. ^ is no longer 
simply £v cos{9 — 4>). First of all, v — v{<j)) and moreover 
V need not be parallel with k. Formally, we could replace 
6* - in Eq. ^ hy 9 - ^(0) with ^ defined by v{k) = 
(cos ^, sin ^)z) ((/)). Single-band equations (|7l8p should be 
replaced by 



cos^(0) = w{(t))a{(j)) 
sin^{(f>) = w{(f>)b{(l)) ■ 



w{cP,<P')a{^') (58) 
«;(0,0')&(0')- (59) 



These two equations for a{(t>) and 6(0) are still completely 
decoupled. With some luck, cos^(0) can be reasonably 
expanded in terms of cosines and sines of and higher 
multiples of but the i'(0) and w(0') terms will most 
likely make an analytical solution of Eq. ([55)) impossible 
for realistic anisotropic Fermi surfaces. The solution is, 
however, not difficult to obtain by numerical means. Af- 
ter discretization of the angular variable into n steps, 
Eq. ()58p constitutes an n x n system of linear equations. 

Once a(0), 5(0) arc known, the non-equilibrium distri- 
bution function is readily written as 

/(fc,£)-/o = /(0,0)-/o= (60) 
-e£w(0)(-9Jo)[a(0)cos6l + h{(j,) sin 6*] . 

Note that the spectral function —defo — S{eF — e^:) de- 
pends now both on k and 0. 

A rather straightforward generalization of Eq. ([50)1 
and the appropriate pair of integral equations ()58l59p to 
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multiband systems is possible. For instance, the anal- 
ogy of the two coupled Eqs. ([50]) for anisotropic band 
structure reads 



cos^+(0) = u'+(0)a+(0)— / d(j)' X 



cosC_(0) = w_ (0) a_ ((/)) 



V -{</>) 



where v± — (cos^±, sin^±)ti±(0) are the Fermi velocities 
of the two bands and the four quantities w±±[4',4>') have 
to be calculated in the spirit of Eq. (I57[) . 
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